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ABSTRACT 

Recently, the Fermi space telescope has discovered two large 7-ray emission regions, the so-called 
"Fermi bubbles" , that extend up to ~ 50° above and below the Galactic center. The 7-ray emission 
from the bubbles are found to follow a hard spectrum with no significant spatial variation in intensity 
and spectral shape. The origin of the emission is still not clearly understood. Suggested explanations 
include injection of cosmic-ray nuclei from the Galactic center by high-speed Galactic winds, electron 
acceleration by multiple shocks and stochastic electron acceleration inside the bubbles. In this letter, 
it is proposed that the 7-rays can be the result of diffusive injection of Galactic cosmic-ray protons 
during their propagation through the Galaxy. Considering that cosmic rays undergo much slower 
diffusion inside the bubbles than in the averaged Galaxy and at the same time suffer from inelastic 
collisions with the bubble plasma, this model can explain the observed intensity profile, the emission 
spectrum and the measured luminosity without invoking any additional particle production processes 
unlike other existing models. 
Subject headings: cosmic rays — diffusion — Galaxy: halo — gamma rays: galaxies 



1. INTRODUCTION 

Recent detailed analysis of the Fermi-LAT data has 
discovered two giant 7-ray emission regions extending 
up to ~ 50° (~ 10 kpc) in Galactic latitude above and 
below t he Galactic ce nter with a width of ^ 40° in lon- 
gitude (jSu et al.ll2010( ). The 7-ray emission regions, now 
popularly known as the "Fermi bubbl es", coincide well 
with the WMAP haze at lo w latitudes (|FinkbeineHl2004l: 
iDobler fc Finkbeinerl l2008f) and share their edges with 
the ROSAT X-ray map (jSnowden et al.iri997f) . More re- 
cently, the Fermi bubbles are also found to be coincident 
with two giant radi o lobes that appear t o originate from 
the Galactic center (jCarretti et al.ll2013[ ). These correla- 
tions found in the multi-wavelength observations seem to 
suggest that the Fermi bubble 7-rays (measured in the 
range of ^ 1 — 100 GeV) are produced by high-energy 
electrons via inverse Compton scattering of ambient low- 
energy photons, as the same electrons can also simul- 
taneously produce radio syi ichrotron radiation s in the 
presence of magnetic fields (jDobler et al.ll2010l ). More- 
over, the fact that the Fermi bubbles are well symmetric 
across the Galactic plane and also well centered on the 
Galactic center encourages to assume the bubbles, and 
also the associated 7-r ay emission, to have their origin 
at the Galactic center (jSu et al.ll2010f ). 

However, because of severe radiative losses, it would 
be extremely difficult to transport high-energy electrons 
from the Galactic center to the far edges of the bub- 
bles. For electrons relevant for producing the Fermi bub- 
ble 7-rays, convection by Galactic wind would require a 
wind speed as high as ~ 10^ km s^^ which is more than 
an order of magnitude larger than the typical Galactic 
wind speed of ~ 200 — 300 km s~^, and diffusive trans- 
port would need a diffusion coefficient which is ~ 2 — 3 
orders of magnitude larger than the standard Galac- 
tic value. This problem can be overcome using models 
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based on jet activity of the central active galactic nucleus 
where j et speeds o f over 10^ km s~^ is rea dily achievable 
(|Guo fc Mathe-^ 120121 : lYang et all 120121 ) . An alterna- 
tive solution would also be to consider production of elec- 
trons inside the bubbles itself, for instance, by generat- 
ing multiple shock waves by perio dic star capture by the 
central supermassive black hole (jCheng et al.l 1201 1[ ) or 
by second-order Fermi acceleration by the plasma wave 
turbulence present throu ghout the volume of the bubbles 
(jMertsch fc SarkaillMl . 

On the other hand, models based on the 7r°-decay ori- 
gin of 7-rays suffer less constraint at least from the par- 
ticle injection point of view. It has been shown that if 
cosmic-ray nuclei from the Galactic center are injected 
into the bubbles by fast wind of speed ~ 10'^ km s~^, 
and if the cosmic rays remain trapped inside the bubbles 
for over ^ 10^° yr, then the inelastic collisions of cosmic 
rays with the low-density bubble plasma, producing tt*^ 
and IT mesons, can expla in both the observed 7- rays 
and the radio emissions (Crocker fc Aharonianll201lD. In 
almost all the models proposed so far (see lSu et al.ll2010l 
for more scenarios), an additional process of high energy 
particle production, either inside or outside the bubble, 
has been considered in order to explain the observed 7- 
rays. In this letter, we present a "minimal" model which 
does not invoke any additional sources or particle pro- 
duction processes other than the sources responsible for 
the production of bulk of the Galactic cosmic rays. 

In our model, it is assumed that cosmic rays (consid- 
ered mainly as protons), after leaving their sources, un- 
dergo diffusive propagation through the Galaxy. If the 
Fermi bubbles are absent of sources, the diffusive stream- 
ing of cosmic rays in the direction of density gradient 
can result into a net flux of cosmic rays injected into 
the bubbles. Inside the bubbles, cosmic rays suffer losses 
due to inelastic collisions with the bubble plasma which, 
in steady state, is assumed to balance the injected flux. 
If the injection is uniform throughout the surface of the 



bubbles, the total amount of cosmic-ray energy injected 
into each bubble per unit time (assuming spherical shape 
based on their projected images shown in lSu et aI]|2010D 
is given by L = AirR^ue, where R represents the radius 
of the bubble, u the streaming or injection velocity and 
e the energy density of cosmic rays. For cosmic rays 
streaming at Alfven velocity va ~ 10^ cm s~^ which is 
the maximum streaming velocity in the case of scatter- 
ing by self-generated Alfven waves, the total cosmic-ray 
power injected is calculated to be L ~ 4 x 10''*' ergs s~^ 
for R = 4.5 kpc and e = 1 eV cm~^, the locally mea- 
sured cosmic-ray energy density. This amount of injected 
power is '^ 2 orders of magnitude larger than the power 
required in cosmic-ray protons to produce the measured 
7-ray luminosity of ^ 2 x 10'^^ erg s~^ from each bubble. 
Even if the cosmic-ray density in the halo is lower by 
one order of magnitude with respect to the local value, 
this rough estimate shows that the injection of some frac- 
tion of Galactic cosmic rays into the bubbles can easily 
account for the measured 7-ray luminosity. 

For cosmic rays with an equilibrium number density Ng 
in the Galaxy and propagating with diffusion coefhcient 
Dg, their injected flux into the bubbles is given by Finj = 
Dg\/Ng ex DgNg, calculatcd at the edges of the bubbles. 
Since Ng is related to the cosmic-ray source spectrum Q 
as Ng ex Q/Dg, we get Finj « Q. Thus, cosmic rays 
streaming into the bubbles will follow the same spectral 
shape as the source cosmic-ray spectrum in the Galaxy. 

Once injected, cosmic rays are assumed to undergo dif- 
fusive propagation inside the bubbles. If the plasma in- 
side the bubbles is extremely turbulent as suggested by 
X-ray observations of the Galactic bulge (jYao fc Wang! 

t07) or if the magnetic field lines are highly tangled 
IcQuinn fc Za ldarriaga"2011"). then the injected cosmic 
rays can take a considerable amount of time to penetrate 
the bubbles due to much slower diffusion. The longer cos- 
mic rays take to penetrate, the more they can be affected 
by inelastic collisions with the bubble plasma resulting 
into depletion of density towards the center of the bub- 
ble. Such a density profile is also suggested by observa- 
tions at low energies. The uniform projected intensity 
distribution of 7-rays observed from the Fermi bubbles 
seem to imply a non-uniform 7-r ay emissivity d istribu- 
tion that peaks towards the edge (|Su et al.ll2010( ) . More- 
over, for a diffusion coefficient that scales with energy, 
the steady state cosmic-ray spectrum is expected to be 
close to Nb (X Finj c>c Q near the boundary, and flatter 
towards the center. This is important because for the 
hadronic origin of Fermi bubble 7-rays, the 7-ray spec- 
trum should mimic the spectrum of the bubble cosmic 
rays. The obser ved 7-r ay spectral index of ^ 2.0 — 2.1 
above ~ 1 GeV (|Su et a l. 2010) is slightly smaller than 
the cosmic-ray source index of ~ 2.1 — 2.3 required to 
produce the locally observed cosmic-ray spectrum. 

2. MODEL CALCULATIONS 

2.1. Cosmic-ray spectrum 

In the absence of sources, the steady state distribution 
of cosmic-ray protons inside a bubble can be described 
by the diffusion-loss equation: 
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where Ni,{r, E) represents the differential number density 
of cosmic rays, E is the kinetic energy, D), represents the 
diffusion coefficient inside the bubble, t{E) = l/{ni,v(j) 
represents the inelastic collision time with the bubble 
plasma of density ri;,, v the velocity of cos mic rays and 
y{E) the inelastic collision cross-section (iKelner et al.l 
|2006[) . In spherical symmetry, Eq. (1) reduces to the 
modified spherical Bessel differential equation as. 
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where we have written x — r/^/rDl with r representing 
the radial variable. The solution of Eq. (2) is readily 
available as. 
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The second term of the solution can be neglected as it 
blows up at the center of the bubble (r = 0), and the 
constant A can be determined using the fiux continuity 
relation at the bubble boundary r = R, 



Db 



dNb 
dr 



= D„ 



dN„ 



r=R, 



dr 



(4) 



r=R 



where the right hand side represents the flux of Galactic 
cosmic rays injected through the surface of the bubble 
F 



in]- The final solution of Eq. (1) is then obtained as, 
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where X = R/^fFlTb. 

To estimate Finj^ we write the Galactic cosmic-ray 
density as a functi on of perpendic ular distance z to the 
Galactic plane as ()Thoudamll2008[ ) . 
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where Q{E) represents the source spectrum, and Dg — 
2.9 X 10^*(i?/i?o)" cm^ s^^ is the diffusion co efficient in 
the G alaxy with £^0 = 3 GeV and a = 0.6 (jThoudaml 
[20M ). The function 
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has a very weak energy dependence, where Ji is the 
Bessel function of order 1, H represents the halo bound- 
ary taken to be large enough to contain the Fermi Bub- 
bles, 3? the radial size of the source distribution taken 
as 20 kpc and 77 the averaged surface density of inter- 
stellar gas in the Galactic plane. Eq. (6) assumes that 
both the sources and the interstellar gas are uniformly 
distributed in the Galactic plane. In terms of the local 
density Ng{E) which is taken to be at z = 0, Eq. (6) 
can be expressed as. 
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is also almost independent of energy. The cosmic-ray 
density gradient along z is then given by, 
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where C = (V/i)//2, hereafter referred to as the injec- 
tion coefficient, will determine the amount of Galactic 
cosmic rays injected into each bubble and depends on 
the distance z. In general, C may vary for different posi- 
tions and also for different directions in the Galaxy. But, 
for the present study, we neglect variations in C and as- 
sume uniform cosmic-ray injection over the surface of the 
bubbles. We then treat C as a parameter which will be 
determined based on the measured 7-ray data. From Eq. 
(10), it can be noticed that as N°{E) oc Q{E)/Dg{E), 
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Eq. (11) shows that the flux of Galactic cosmic rays en- 
tering the bubbles Finj — Dg\7Ng j^^^^ follows the source 
cosmic-ray spectrum in the Galaxy as mentioned before. 
We further assume that the diffusion coefficient inside 
the bubble scales with the galactic value as Dj, — KDg, 
where if <C 1 is a constant. And, the cosmic-ray source 
index is taken to be F = 2.2 so that, for the assumed Dg, 
the demodulated local cosmic-ray intensity {v/Ati)N^ re- 
produces the measured proton spectrum between ^ 10 
GeV and 1000 GeV (Adriani et ah 2011). 

2.2. Gamma-ray em,ission 

In the present model, 7-rays from the bubbles are con- 
sidered as the decay product of 7r° mesons produced 
from the inelastic collision of cosmic rays with the bub- 
ble plasma. For a constant plasma density inside the 
bubbles, the 7-ray emissivity is expected to follow the 
distribution of cosmic rays. The 7-ray emissivity as a 
function of the bubble radius is calculated using 
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where E^ denotes the 7-ray energy, E^o and 1 
note the total and the rest mass energy of the pion re- 
spectively, and i?™" = Ery + mlo/AE^. The tt" mesons 
emissivity is given by (jKelner et al.ll2006l ) 

q^o{r,E^o) = jcnbaiT)Ni,{r,T) (13) 

k 

where c is the velocity of light, and Ni,{r, T) is the radial 
dependent cosmic-ray proton density given by Eq. (5) 
but taken as a function of the total energy T which is 
related to the pion energy as T = mp + E^o/k. Following 
iKelner et al"] (|2006( ). we take k = 0.17, n == 1, and 
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where L = ln(r/l TeV) and Tth = 1.22 GeV is the 
threshold energy for 7r° meson production. 
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Fig. 1. — Normalized distribution of 10 GeV cosmic rays inside a 
Fermi bubble for different values of D5 with K = (0.2 — 5) X lO""^ 
and rui = 0.01 cm""^. The thin solid line is for 1 GeV particles for 
the case of X = 6 x 10"''. 

The 7-ray intensity in a given direction characterized 
by the Galactic longitude I and the latitude b is calcu- 
lated as 
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where the integration is performed along the line of sight 
distance y, and the integration limits are determined 
from the points of intersection of the line of sight with 
the bubble surface. In calculating Eq. (15), the emissiv- 
ity previously calculated as function of the bubble radius 
r has been carefully expressed as function of y. 

3. RESULTS AND DISCUSSIONS 

The distribution of cosmic rays inside the bubbles is 
governed by the competition between diffusion and in- 
elastic collision. Diffusion tends to uniform the distribu- 
tion while inelastic collision does the opposite by remov- 
ing particles. A slow diffusion process will be more af- 
fected by inelastic collisions, leading to more exclusion of 
cosmic rays from the bubble interiors. Thus, faster diffu- 
sion will produce more uniform distribution than slower 
diffusion. 

For energies of our interest, the inelastic collision cross- 
section is known quite accurately and hence, for a given 
plasma density, the inelastic collision time t{E) can be 
fixed. On the other hand, the diffusion coefficient inside 
the bubbles is largely unknown although a value smaller 
than the Galactic averaged is suggested (|Yao fc Wand 
|2007() . Figure 1 shows the normalized distribution of 10 
GeV cosmic rays (thick lines) inside a bubble for different 
values of Db with K = (0.2 - 5) x 10"^ taking nt = 0.01 
cm~^ ()Su et al.ll2010D . It can be seen that the distribu- 
tion becomes flatter towards the center of the bubble as 
K takes larger values which is due to the shorter diffusion 
times t ^ R^ /Db of cosmic rays to reach the center. Also 
shown for comparison is the distribution of 1 GeV parti- 
cles (thin solid line) for the case oi K = 6 x 10""*. The 
distribution is steeper than the corresponding distribu- 
tion of 10 GeV particles (thick solid line). In our model, 
the cosmic-ray distribution is expected to be more uni- 
form at higher energies due to the energy scaling of the 
diffusion coefficient as Db oc E^-^. 
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Fig. 2. — Projected 7-ray intensity profile of the Southern bubble 
for the averaged 1 — 5 GeV (solid line) and 5 — 20 GeV (dashed 
line). Also shown are the predictions for 100 GeV (dotted line) and 
1000 G eV (dot-dashed line) energies. Data are taken from lSu et al.l 
l(20TOl 'l. 

Once the value of r is fixed, the cosmic-ray distri- 
bution inside the bubbles is determined by the choice 
of Db (through the constant parameter K). Here, its 
value is chosen such that the resulting projected 7-ray 
intensity distribution matches the measured profile. It 
is found that choosing K = 2.2 x 10"'' produces a good 
fit to the measured data as shown in Figure 2, where 
the data represents the measurements from the South- 
ern bubble. The calculation assumes R = 4.5 kpc and 
the bubble center t o be at 5 kpc above the Galactic cen- 
ter ()Su et al.ll201Clf) . The model predictions for both the 
averaged 1 — 5 GeV (solid line) and 5 — 20 GeV (dashed 
line) are added with backgrounds obtained by fitting hor- 
izontal lines to the respective data between 5° and 40°. 
It can be mentioned that neither the ha dronic model pre- 
sented in iCrocker fc AharonianI ()2011l) nor the leptonic 
mode l based on diffusive shock acceleration ([Cheng et al.l 
I2OIID can satisfactorily explain the measured sharp edges 
shown in Figure 2. Both these models predicted a con- 
stant volume emissivity throughout the bubbles which 
will produce soft edges on the projected profile. Also 
shown in Figure 2 are the predictions for 100 GeV (dotted 
line) and 1000 GeV (dot-dashed line) 7-rays which can be 
tested in future. Their profiles are relatively steeper to- 
wards the edge than at low energies, and also show softer 
edges, reflecting the more uniform particle distribution at 
higher energies. It can be noted that our result at high 
energies, say at 1000 GeV, is clearly different from that 
expected fro m the leptonic stochastic acceleration model 
presented in lMertsch fc Sarkaii ()201lD which predicted a 
significant edge brightening at high energies. This dif- 
ference in the high energy radial profile can be used to 
differentiate between the two models in future. 

The 7-ray spectrum averaged over a whole bubble is 
shown in Figure 3 (top). The calculation assumes F = 2.2 
as mentioned before and an injection coefhcient of C = 
4.6 X 10~^ kpc~^. The latter corresponds to an injection 
flux of 0.46% of the local cosmic-ray density injected with 
velocities of ~ 10^ cm s~^ at few GeVs into each bubble. 
Note that this is also the same injection strength used 
to produce the results presented in Figure 2. The model 
prediction is found to be in good agreement with the 
data in the energy range of ~ 1 — 100 GeV where the 
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Fig. 3. — Gamma-ray spectra for a whole bubble (top), and for 
the bubble inner and outer regions (bottom). T he outer r e gion i s 
taken as a shell of 1 kpc thickness. Data are from lSu et al.l II2010I V 

measurement uncertainties are small. It is interesting to 
see that the same source index required to explain the 
measured cosmic-ray spectrum also reproduces the Fermi 
bubble 7-ray spectrum. 

Spectra for t he inner a nd ou ter regions of a bubble 
separated as in iSu et al.l ()2010l ) are shown in Figure 3 
(bottom). The outer region is taken as a shell with thick- 
ness 1 kpc. The present model produces a steeper 7-ray 
spectrum in the outer region than in the inner region. 
For the inner, the model prediction agree nicely with 
the measured data between ~ 1 — 100 GeV, but for the 
outer region, the model spectrum seems to be somewhat 
steeper compared to the data. Choosing a smaller diffu- 
sion index inside the bubble will produce more spectral 
uniformity in the two regions. However, for any reason- 
able value of the diffusion index, our model prediction is 
expected to be different from that of the leptonic model 
(JMertsch fc Sarkadl2011| ) in which a steeper spectrum in 
the inner region is expected. 

In our calculations so far, we have assumed that the 
cosmic-ray diffusion index inside the bubbles is the same 
as that in the Galaxy. It is possible that the index is dif- 
ferent in the two regions depending on the nature of the 
turbulent wave responsible for scattering the cosmic-rays. 
Actually, a higher turbulent region is expected to have a 
steeper turbulent wave spectrum w{k) ex fc~" (where k 
represents the wave number) which then translates into 
a flatter diffusion index. For a Kolmogorov type wave 
spectrum which is characterized by a turbulent index of 




Data: Averaged (1 -2)&(2-5) GeV n 
Data: Averaged (5-1 0)&(1 0-20) GeV o 

Model: Averaged (1-5) GeV 

Model: Averaged (5-20) GeV 

Model: 100 GeV 

Model: 1000 GeV 



V^*^*^^'^**'^ 



10 20 

Degrees from the edge 



Data: FERMI 0.5-1 GeV IC template (Inner bubble) • 
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Fig. 4. — Same as Figure 2 and Figure 3 (bottom panel), but for 
Di, (X £1/3 and K = 5.5 X 10""'. 



a — 5/3, the diffusion index becomes a = 1/3. We have 
also carried out calculations taking a = 1/3 inside the 
bubbles (keeping a — 0.6 for the region outside). For the 
same value of injection coefficient obtained above for the 
case of a = 0.6, it is found that a value oi K = 5.5 x 10""* 



is required in order to explain both the measured inten- 
sity radial profile and the measured spectrum simultane- 
ously. The results are shown in Figure 4. For the radial 
profile, noticeable changes are observed only at high en- 
ergies. For instance, the profile for 1000 GeV becomes 
comparatively fiatter with shaper edges than the corre- 
sponding result shown in Figure 2 for the case of a = 0.6. 
At low energies, the radial profiles look very much simi- 
lar in the two cases. Also, from Figure 4 bottom panel, it 
can be noticed that the difference in the predicted spec- 
tra between the inner and outer regions of the bubble 
becomes less. Thus, more sensitive spectral measure- 
ments at high energies can give a better constrain on the 
diffusion index in the present model. 

It should be mentioned that although our simple model 
can successfully explain the observed radial profiles for 
the Southern bubble, it cannot explain the dip at ^ —12° 
observed in the 5 — 20 GeV profile of the Northern bubble 
()Su et al.ll2010[ ). Moreover, also the apparent hardening 
of the observed spectrum below ~ 1 GeV cannot be satis- 
factorily explained (see Figures 3 and 4 bottom). So far, 
none of the existing models can explain this spectral be- 
havior. Future accurate measurement of this low-energy 
spectral shape will be crucial for better understanding 
the nature of high-energy particles responsible for the 
7-ray emission. 

4. CONCLUSIONS 

In this letter, we have proposed a minimal model for 
the Fermi bubbles which, unlike other existing models, 
does not consider any additional particle production pro- 
cesses or sources other than those responsible for the pro- 
duction of the Galactic cosmic rays. The model consid- 
ers that the 7-rays from the bubbles can be the result 
of injection of Galactic cosmic-ray protons during their 
diffusive propagation through the Galaxy. It is found 
that this simple model can explain many of the observed 
properties of the Fermi bubbles. 

I wish to thank D. Jones and A. Achterberg for insight- 
ful discussions. I also thank J. R. Horandel, H. Falcke, 
and P. Groot for constant encouragement. 
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